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Abstract 

For each even N > Cdt d we prove the existence of a well separated 
spherical f-design in the sphere S d consisting of N points, where Cd 
is a constant depending only on d. 
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1 Introduction 

In this paper we will discuss the interrelation between several classical opti- 
mization problems on spheres S d such as minimal equal-weight quadratures 
(spherical designs), best packing problems, minimal energy problems. Typi- 
cally for d = 1 the optimal configuration for all of these problems is a regular 
polygon. However for d > 2 exact solutions are known in very few cases. 
Even asymptotically optimal solutions sometimes are very hard to obtain, 
see for example Smale's Problem 7 for 21st century |18j . 

We will prove the existence of certain configurations in S , which are 
spherical t-designs with asymptotically minimal number of points, and si- 
multaneously have asymptotically the best separation property. These con- 
figurations also provide approximate solutions for several other optimization 
problems. 
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Let S d = {16 ]R d+1 : |x| = l} be the unit sphere in M. d+1 equipped 
with the Lebesgue measure fid normalized by fid(S d ) = 1. A set of points 
xi, . . . , %n € S d is called a spherical t-design if 



for all polynomials in d+1 variables, of total degree at most t. The concept of 
a spherical design was introduced by Delsarte, Goethals, and Seidel [9]. For 
each d, t G N denote by N(d, t) the minimal number of points in a spherical 
t-design in S d . The following lower bound 



is proved in [9], see also the classical monograph [8]. Let us give a brief 
history of asymptotic upper bounds on N(d, t) for fixed d and t — > 00. First, 
Seymour and Zaslavsky [17] have proved that spherical designs exist for all 
positive integers d and t. Then, Wagner [19] and Bajnok [3] independently 
proved that N(d,t) < C d t Cd4 and N(d,t) < C d t Cd \ respectively. Korevaar 
and Meyers [12] have improved these inequalities by showing that N(d, t) < 
Cdt (S+d)/2_ They have algo con j ectured that N(d,t) < C d t d . Note that © 

implies N(d,t) > c d t d . In what follows we will use the notations C d , L d , etc. 
(c d , \d, etc.) for sufficiently large (small) constants depending only on d. 

Korevaar and Meyers were motivated by the following problem coming 
from potential theory: How to choose N equally charged points x%, ...,x^ in 
S 2 to minimize the value 



The classical Faraday cage phenomenon states that any stable charge dis- 
tribution on the compact closed surface cancels the electric field inside the 
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(1) 





r G (0,1)? 



surface. According to this model the minimal value of U r should rapidly 
decay to 0, when N grows. 

It was shown in (12] that if the set of points x±,...,xn is a spherical 
i-design, where t > cN 1 / 2 then 

U r (x u ...,x N )<r c > Nl/2 . 

The estimate is exact up to the constant in the power. 

In we have suggested a nonconstructive approach to obtain an opti- 

mal asymptotic bound for N(d, t) based on the application of the topological 
degree theory. We have proved the following 

Theorem A. For each N > Cdt d there exists a spherical t-design in S 
consisting of N points. 

This proves the Korevaar-Meyers conjecture. 

Now we will give the definition of a well separated sequence of configura- 
tions. The sequence of A-points configurations = {xm, • xnn} in S d 
is called well separated if 

(2) min \x iN - x jN \ > A d A~ 1/fl! , 

for some constant A^ and all N > 2. The inequality fl2]) is optimal up to the 
constant A^. That is, there exists a constant such that for any Appoint 
configuration {xi, ...,xn} in S d there exists i ^ j such that 

\xi-Xj\ <L d N- l ' d . 

Many authors have predicted the existence of well separated spherical t- 
designs in S d of asymptotically minimal cardinality 0(t d ) as t — > oo, see, 
e.g. [2] and [TTJ. Moreover, in pTJ was shown that if such spherical designs 
exists then they have asymptotically minimal Riesz s-energy. In this paper 
we prove the existence of above mentioned spherical designs. Our main result 
is: 

Theorem 1. For each d > 2 there exist positive constants Cd and A^ de- 
pending only on d such that for each t £ N and each even N > Cdt d 



3 



there exists a spherical t-design in S consisting of N points {xi]f =1 with 



3lj it 3u i 



>\ d N- 1 l d fori^j. 



Theorem [I] is a natural generalization of Theorem A. The paper is orga- 
nized as follows: In Section [2] we will reduce Theorem [1] to the construction 
of a certain iV-tuple of maps x 1; . . . , x^r : Vt — > S d . Then in Section 3 we will 
prove several auxiliary results concerning area regular partitions of sphere S d 
and spherical Marcinkiewicz-Zygmund type inequalities. Finally in Section 4 
we will construct the maps Xj, . . . , Xjv, proving Theorem 1. 

2 Application of topological degree theory 

Let Vt be the Hilbert space of polynomials P of degree at most t on S d such 
that 

P(x)dfi d (x) = 0, 



equipped with the usual inner product 



(P,Q)= / P(x)Q(x)dfi d (x). 
Js d 

By Riesz representation theorem, for each point x e S d there exists a unique 
polynomial G x G Vt such that 

(G x ,Q) = Q(x) for all Q E V t . 

Then a set of points Xi, . . . ,Xn £ S d forms a spherical t-design if and only if 

(3) G X1 H h G XN = 0. 

The gradient of a different iable function / : — y R is defined as 
df fdf df \ 

z = (£i,...,Un). 



dx ' \d£i' ' d£ d+1 
For a polynomial Q EVt we define the spherical gradient 

d „ / x 



v «w : =s«(e)- 
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Now we will use the following result from the topological degree theory, see, 
e.g., ! 1 •") . Theorem 1.2.6 and Theorem 1.2.9]. 

Theorem B. Let f : MJ 1 — > W 1 be a continuous map and Q an open bounded 
subset, with boundary dQ, such that G fl C M. n . If (x,f(x)) > for all 
x G dQ, then there exists i6fl satisfying f(x) = 0. 

We will apply Theorem B to the following open subset of a vector space Vt 
(4) n= \ pe V, 



\VP(x)\d/i d (x) < 1 

s d 

Observe that if continuous maps Xj : V t ->■ S d , i = 1, . . . , iV satisfy for all 

P edn 

N 

^P( Xl (P))>0 

i=l 

then there exists a spherical t-design on S d consisting of N points. To this 
end let us consider a map / :V t — » Vt defined by 

f(P) = G Xl (p) H h G XN ( P ). 

Clearly 

(p,/(p)) = f;p(x i (p)) I 

i=l 

for each P G Vt- Thus, applying Theorem B for the map /, vector space 
Vt, and the subset fl defined in PJ immediately gives us the existence of a 
polynomial P G Q such that /(P) = 0. Hence, by (j3J), the images of this 
polynomial xi(P), ...,x/y(P) form a spherical t-design in S" d consisting of 
points. If additionally there exists a constant such that |xj(P) — Xj(P)| > 
ArfiV _1//d for all % j, and Fefi, then the above mentioned spherical t-design 
will be well separated, proving Theorem [TJ 

The maps Xj, i — 1, . . . , N, will be constructed in Section H] below. 
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3 Auxiliary results 



For x,y G S d denote by dist(x,y) = arccos((x, y)) the geodesic distance 
between x and y. Also for the set A C S define the geodesic distance 
between x and A as follows 

distfx, A) = inf dist(x, y). 

Recall that a spherical cap of radius r with center at x G S is the set 
A(x, r) = {z G S d \ dist(x, z) < r}. 

Below we will extensively use the notion of an area-regular partition. Here 
is a definition. 

Let 7Z = {Ri, . . . , Rn} be a finite collection of closed sets Ri C S d such 
that Lif =1 Ri = S d and fi d (Ri n Rj) = for all 1 < z < j < iV. The partition 
7?. is called area-regular if (j>d{Ri) = i = 1, ■ ■ ■ ,N. The partition norm 
for 1Z is defined by 

1 172.1 1 = maxdiami? = max max distfx, y). 

R&TZ ReTZ x,y£R 

It is easy to prove using isodiametric inequality that each Ri has diameter at 
least CdN~ l l d . Therefore, ||72|| > CdN^ 1 ^ 1 for each 1Z. However for some d 
and each iV G N there exists an area regular partitions 71 = {R\, . . . , Rn} of 
S d with diameter at most CdN~ l / d . Such area-regular partitions are used for 
many optimization problems where a well distributed set of iV points on a 
sphere having no concentration points is needed, see e.g. [TJ, [B], [2], and [T5] . 
In this paper we need area regular partitions of small diameter with additional 
constraint of geodesic convexity. A subset A C S is geodesically convex if any 
two points x, y G A can be connected by a geodesic arc contained in A. The 
partition 71 = {Ri, Rn} is said to be convex if each set R i: i = 1, . . . , N, 
is geodesically convex. First will prove the existence of convex area regular 
partitions of "small" diameter. 

Theorem 2. For each even N G N there exists a convex area-regular parti- 
tion TZ = {Ri, . . . , Rn} such that \\7Z\\ < K d N~ l / d , where K d is a constant 
depending only on d. 
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The following construction for the sphere S 2 and N = 6n 2 , where n £ N, 
is given by Alexander in pQ . Let us first explain his simple and elegant proof. 

We begin with a spherical cube, and consider one of its faces. Using n — 1 
great circles from the pencil determined by two opposite edges we can cut the 
face into n slices of equal area. Each slice can be cut into n quadrilaterals of 
equal area using great circles in the pencil in the other pair of opposite edges 
of the face. The diameters of the quadrilaterals are of the right magnitude. 

This construction has an obvious generalization to higher dimensions. 
Start with the appropriate spherical hypercube, then divide each face into n 
equal pieces, and so on. In this way we obtain a partition of S into 2(d+l)n d 
regions. 

Now our goal is to generalize the approach of Alexander and prove The- 
orem [2] for all even N £ N. For each m £ N and a vector with positive 
coordinates a = (a 1; . . . , a m ) £ denote by P(a) the m-dimensional rect- 
angle [— ai,ai] x ... x [—a m ,a m ]. Also, a measure drj(x) = a(x)dx defined 
on P(a) is said to be M-uniform if 



and drj be a M-uniform measure defined on P(a). Then for each N £ N 
there exists a partition of P(a) into N m-rectangles Pi, ... , Pn with facets 
parallel to the corresponding facets of P(a) such that 



a(x) < Ma(y) for all x, y £ P(a). 



To prove Theorem |2] we need the following lemma. 



Lemma 1. Let a £ MJ? with 



max — < B 

l<i,j<m aj 





and 



(6) 



diamPi < C(m,B,M)a 1 N~ 1/ 



for alii — 1, . . . , 



N. 
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Proof. We will prove the lemma by induction on m. For m = 1 first we will 
choose a point t\ G [ — 0\ , CLi 1 such that 

M x ) = T7 / d V( x )- 
Then similarly we choose i 2 G °i] such that 

and so on. Finally, we get the partition of P(a) = [— ai, ai] into A" segments 
Pi = [— ai, ti], P2 = £2],- • •; Pv = [^JV-ij a i] satisfying (jSj) by its definition. 
Moreover, M-uniformity of 77 implies ([6]) with C(l, B, M) = 2M. Assume 
that the lemma is true for m = I — I. Let us prove it for m = I. Put k = [N 1 ' 1 ], 
s = [N/k], and r = N mod k. Also for a fixed a = (a 1; . . . , ai) G denote 
by Pi (a) the rectangle [— ai,aj x . . . x [— a/_i, aj_i]. To obtain a necessary 
partition of P(a) first we choose step by step points — a n = t Q < t\ < . . . < 
tk-i < tk = a n such that 

J drj{x) = ^— / dq(x) 

J[ti-T_M]y-Pl(.a) iV JP(a) 

for all i — 1, . . . , r, and 

[(.-iA]xPi(a) iV ip(a) 

for alH = r + 1, . . . , k. Clearly, 

20 + 1) ,/, 

(7) \U - < 1 - > Ma l < AN I BMa\. 

Consider the following measures on Pi (a) 

drji(x) = a>i(x)dx, 

where 

(8) ati(x) = a i (x 1 , . . . ,2^1) = / a(x 1 , . . . , xi)dx t 
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for all (x\, . . . , Xi-i) G -Pi (a) and % = 1, . . . , k. Clearly each rji is 2M 2 -uniform. 
Hence, by induction assumption for each i = 1, . . . k and iVj e N there exists 
a partition of Pi into N{ rectangles P^i, . . . , Pi,jv, with facets parallel to the 
corresponding facets of Pi (a) such that 

(9) / dr)i(x) = — I drn(x), 

JP l}j iV i JP^a) 

and 

(10) diamPij < C(l - 1, B, M) ai NJ l/{l ~ l) 

for all j = 1, . . . , iVj . Choose iVj = s + 1 for i = 1 , . . . , r and Ni = s for 
i = r + 1, . . . , k. Consider the following partition of P(a) into N rectangles 

Pij x [U-uU], i — I, ■ ■ ■ ,k, j — 1, . . . 

By (jSJ) and (jUJ) we immediately get that 

rfr/(x) = / dr)i(x) = -rr dr)i{x) = -J- / dr/(x) 

for alH = 1, . . . , k and j = 1, . . . iVj. So, for this partition (jSJ) holds. Finally, 
combining ([7]) with (jTt)]) we get (J6]) for some constant C(l, B, M). Lemma [1] 
is proved. □ 

Now we are ready to construct a required convex area regular partition. 

Proof of Theorem^ We may assume that N > 8c? 2 . For each d+ 1-rectangle 
P(a), a G S d , denote by i^i-i its facet X{ = a>i and by F21 its facet Xi = — dj, 
i = 1, . . . , d + 1. One can naturally associate with P(a) a convex partition 
{-Ri(a), . . . , R2d+2( a )} of S d , where -Rj(a) = g(Fi(a)), and = for 
all x G \ {0}. 

Consider a one-parametric family of d + 1-rectangles P(a\), where 



aA= ( A 'VV ! '---'V L ^-) 
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Now we will choose such a A = A (AT) that our required convex area regular 
partition could be obtained as a subpartition of {Ri(a x ), . . . , R2d+2( a x)}- 
Consider the function G(X) = (j, d (Ri(a\)) . Clearly, G(l/y/d+ 1) = l/(2d + 
2) (in this case P(a\) is cube). On the other hand we can estimate G\ 
from above by the area of corresponding spherical cap. We have G(X) < 
Hd(A((l, 0, . . . , 0), arccos A)). By an easy calculation we get that 

, m N/2-d([N/(2d + 2)] + l) 
/^(A((l,0,...,0),arccos(l - l/10d)) < — 

for all d G N and N >8d 2 . We see that for N > 8d 2 
N/2-d([N/(2d+2)]+l) 1 



N 



N ~ 2d +2 

Hence, by continuity of G there exists 

(11) A e [1/Vd+ 1,1 - l/10d] 

such that 

_ iV/2-d([iV/(2d + 2)] + l) 

]v ■ 

By symmetry arguments fJ,d(Ri{a\)) = /^(-^(^a)) and 

r ^ n? r u [iV/(2d + 2)] + i 

For each i = 1, . . . ,2d + 2 consider the unique continuous measure rji on 
Fj(aA) such that rji(E) = jjLd(g(E)) for each measurable set E C Fj. Clearly, 
each ?7i is M^-uniform for large enough M^. Choose 

Jvi = ^-^AM + 2)] + l) fa . = 1[2 

iy,= ^' M ; 2 " + 1 fo, < = 3,..., M + 2. 

Now applying Lemma [1] for each <i-rectangle Fi(a\) with measure rji, i = 
1, . . . , 2d + 2 we can get corresponding partition of Fi(a\) into Ni rectangle 
Pij such that 
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for j = 1, . . . , iVj. Moreover 

(12) diam Pij <C d N- 1/d . 

By its definition fi d (g(P id )) = 1/N, i = 1, . . . , 2d + 2, j = 1, , iV<. Now 

we observe that each g(P{j) is a geodesically convex closed set. Indeed, the 
image under the map g of a line segment contained in Fi ia a geodesic arc on 
sphere S d . Therefore, the image of the convex set P^j is geodesically convex. 
Finally, the estimate diam g{Pij) < CaN- 1 ^ follows from (TU]) and §F2§. □ 

Remark. The fact that TZ is convex easily implies that each Ri, i = 1, .., N, 
contains a spherical cap of radius bdN~~ l l d . 

The following Theorem C states that an arbitrary large enough and well 
distributed set of points is "almost" equal weight quadrature formula in S d , 
see [TU Theorem 3.1]. 

Theorem C. There exist constants r<i and Nd such that for each integer 
m > Nd, each rj G (0, 1), an arbitrary convex area-regular partition 1Z = 
{Ri,...,Rn} with \\7Z\\ < rj—, and each collection of points Xi G Ri, i = 
1, . . . , N, the following inequalities 

r l N r 

(13) (1-17)/ \P(x)\df, d (x)<-J2\ P (^)\<a + ri) / \P(x)\d^ d (x), 

Js d ^ ~l Js d 

hold for all polynomials P of total degree at most m. 
To prove Theorem [T] we need the following lemma. 

Lemma 2. For an arbitrary convex area-regular partition 1Z = {Ri, . . . , Rn} 
with \\1Z\\ < ?7^fj, any two collections of points Xi, G Ri, i = 1,...,N, 
the following inequalities 

i - r 

( 14 ) Tt E I VP ^) " VP ^\ < 4 ^ / I VP(x)\dfi d (x), 
N l^i Js* 



(15) 



N 

l-4d7j) / I VP{x)\ dfx d (x) < i- V I VP{xi)\ < (l+Adrj) I | VP{x)\ dfi d i 
Js d ^ ~l Js d 



x), 
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hold for all polynomials P of total degree m > N d . The constants r^ and 
are given by Theorem C. 



Proof. First we will prove (JHJ). Since |VP| = y Pf + . . . + P% +1 , where 
Pi G Vm+i for i = 1, . . . , d + 1, we have 

1 N 1 N d+l 

( 16 ) iv E > VP ^) - VP ^)i ^ iv E E i p i(^) - p i(w)i- 

i=l i=l j=l 

Now, for each Q G ? m+ i we will estimate the value 

i N 

-J2\Q(xi) - Q( yi )\. 
i=i 

Let Ji be the set of indexes i = 1, . . . ,N such that the value Q(x) has the 
same sign for all x G Ri, and /2 be the set of all other indexes, that is the set 
of i — 1, . . . , N, for which there exists a point x G Ri with Q(x) = 0. Let I3 
be the set of indexes i — 1, . . . , N such that > \Q(yi)\. Put Z{ := x iy 

if i G I3 and Zj := ?/i otherwise. Put ^ := if i G ii fl I3, and tj := £j, if 
i & h \ I3. For i G ^2, let t» be a point in i?j such that Q(ti) = 0. We have 

1 N 1 1 1 

< ^ E - ^ E + 1 E 1^)1 - 1 E 1 w)i 

i€/i iG/i ieh i£h 

2 N 2 N 

^ n E \qm\ -^E iwi, 

i=l i=l 

where Zj, G -Rj, i = 1, . . . , N. Thus, by ( 1T3|) we have 

_V|g(x i )-g(y i )|<2r 7 / \Q(x)\dfi d (x). 
So, the inequality ( !T6l) implies 

— J2\VP(xi)-VP(yi)\ <2r]J2 \ V P^d^x) < Ad V / \VP(x)\dfi d (x). 
N 1^1 l^i J s d J s d 
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This proves (I14p . Now by the mean value theorem there exist yi G Ri such 
that 

1|VP(2/0|= / IVP^)!^^), i = l,...,JV. 

J Ri 

Finally, the inequality ( |T5|) follows from f fl^l) and the following easy inequal- 
ities 

N N N 

^ E |VPfa)| - 1 E |VP(x 4 ) - VP( yi )| |VP(x,)| 

i=l i=l i=l 

JV JV 

1=1 1=1 

□ 

The following lemma is crucial to construct of the maps xi, . . . , x./v : Vt — > 
S d in the next section. 

Lemma 3. For x G S d denote by T x the space of all vectors y G M. d+1 with 
(x, y) = 0. Let R G S d be a closed geodesically convex set with diamP < 7r/2. 
Then for each interior point x G R and y G T x \ {0} the following holds: 

(i) there exists a unique x max G R with (x max ,2/) = max 2eJ j(z, y); 

(ii) the map M x : T x \ {0} — » i? jiuen by y x max zs continuous on T x \ {0}; 
fmj /or eac/i w E R and a geodesic 7 : [0, 1] — >• i? wtt 7(0) = x max; 7(1) = w 
the function (y,j(h)) is decreasing on [0,1]. 

Proof. Consider an orthogonal projection p : M. d+l — > T x given by 

p(z) — z — (x, z)x. 

Clearly, 

(17) (z,y) = (p(z),y), 

for all z G R. 
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Denote by S = p(R) the image of R under the projection p. Since 
dist(x, z) < 7r/2, z G R, then p is a homeomorphism between R and S 1 
and the inverse map is given by 

p~ l {u) = u + a/1 — |w| 2 x, u E S. 

Now we will show that 5 is a strictly convex subset of T x , i. e. for any u, 
v £ S, where u ^ v , and /i G (0, 1) the point hu + (1 — /i)t> is an interior 
point of S 1 . To this end we note that 

(18) p~ l {hu + (1 - h)v) = hp" 1 ^) + (1 - h)p~\v) + ax, 

where a > 0. We will use the following simple statement: 
If Wi, W2 G 5 ,d such that (wi, u^) > 0, and W3 = a.\Wi + ^2^2 G S* d for some 
ai,a2 > 0, then w 3 lies on the shortest geodesic connected W\ and wi- 
This statement and the fact that G R immediately imply 

that z/\z\ G R, where z = hp~ l {u) + (1 — h)p~ l {v). Hence, applying again 
the statement for z/\z\ and x we get by (fT8|) that p~ x {hu + (1 — h)v) G R. 
Moreover, since x is an interior point of R, and a > 0, then p~ 1 (hu+ (1 — /z)f ) 
is an interior point of R as well, and therefore hu + (1 — is an interior 
point of 5. 

To prove (i) we will use the known fact that a nonconstant linear function 
given on a closed strictly convex subset in M. d attains its maximum in the 
unique point. Using this fact we get that there exists a unique z m£LX G S such 
that (z max ,y) = m&x zeS (z,y). Finally, by <$TQ we get that a; max = p _1 (z max ). 

Now we will prove (ii). Since p is a homomorphism it suffices to show 
that the composition map p o M x : T x \ {0} — > S is continuous. Note that 
(y,p oM x (y)) = m&x zeS (y, z). Since S is a closed strictly convex set then for 
each e > there exists 5 = 5(e) > such that for all v G T x with \v — y\ < 5 
the diameter of the set {z G S\(v,z) > (v,p(M x (y)))} is less than e. Hence, 
\p(M x (y)) —p(M x (v))\ < e. Thus, the map po M x is continuous at y, and so 
is M x . This proves (ii). 

Finally, we prove part (iii) of the lemma. Let G be a great circle con- 
taining x max and w. There is a unique point w max G G such that (y, u> max ) = 
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max z( zG{y, z )- Now for each z G G we have 
(19) (y,z) = (y,w max )(z 

Hence the scalar product (y, z) is increasing on the both geodesic arcs con- 
necting — U7 max and w max . The geodesic 7 : [0, 1] — > R is an arc of G. To 
prove (hi) it is enough to show that both —w max and w max are outside of the 
arc 7. The point w max is outside of 7 by the definition of x max . Moreover 
(y, imax) > (l/,^) = 0. Therefore, substituting z = x max to ( |T9l) we see that 
(a; ma x,w max ) > 0. Hence, 

dist(x max , -w max ) > tt/2. 

Finally, the fact that diamP < n/2 implies that — w max is outside of 7 as 
well. Thus, the function (y,j(h)) is decreasing on [0, 1]. □ 

4 Proof of Theorem 1 

In Section 2 we explained that it is enough to construct an iV-tuple of con- 
tinuous maps x.i, . . . , Xjv : Vt — > S d such that 

1 N 

-£P(x 4 (P))>0 

for all P G dVt and 

dist (xi(P), x,(P)) > AdiV- 1 ^, 1 < i < j < N, 

for all Pffi, where Q is given by (HI). 
Fix e,5,T] > 0. Consider the function 

{*/e if * < e, 
1 otherwise. 

Let 1Z = {Pi, . . . , Rn} be an area regular partition provided by Theorem [2] 
with N > t d (^j) d . For each i — 1, . . . , N choose a point X{ G Pj such that 
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Ri contains a spherical cap of radius b^N l l d with center at x^. Recall that 
II-RII < KaN- 1 ^. 

Let P E Vt- By Lemma |3] for each i = 1, . . . , N there exists a unique 
Zi = Zi(P) E Ri satisfying 

(20) VP(xi)) = max(x, VP(z<)), 

provided that VP(xj) 7^ 0. In the case VP(xj) = put ^ = Xj. Let 
7[xi,2i] : [0, 1] — > Ri be a geodesic connecting Xj and 2*. Note that I7L. z .i(h)\ = 
dist(xi, Zi) for ft, G (0, 1). Define 

(21) x l (P):= 7 [^]((l-5)^(|VP(x i )|)). 

By the definition of g £ the map Xj, z = 1, . . . , N, is continuous in all P such 
that VP(xj) = 0. On the other hand, if VP(xj) 7^ then Xj is continuous 
in P by Lemma [3] (ii). Thus the maps Xj., . . . x^r are continuous in Vt- The 
following Lemma H] will finish the proof of Theorem [TJ 

Lemma 4. There exist constants e, 5, rj depending only on d such that for 
N > t d (^) d the N -tuple of maps x 1 ,...,x N : V t S d defined by (|2I|) 
satisfies the following properties: 

1 N 

(22) ^E P ( X ^ P )) >0 
/or a// P G <9f2 and 

(23) dist (xi(P), x,(P)) > A d iV- 1/d , l<i<j<N, 
for all PeVl. 

Proof. Fix P E Vt. For each z = 1, . . . , iV choose z^ max E Ri such that 

P(z i,max) — m & x xG-Ri P(x). Denote : — 7[a;i,z<] 

(o e (|VP(xi)|)), where Zj and 
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7 are as in (12"T|) . We can split the sum (122]) into four pieces 

i N ^ N 



i=l i=l 



1 - 

]T(P(^) - P(^, max )) 



N 

i=l 



1 * 



i=l 
i=l 

We will estimate each of these sums separately. 
Clearly, 

N N 

( 25 ) /v 5Z P ( 2 ™) = Z) / ( P (^.— ) - P(x))<W*)- 

1=1 1=1 

Now note that if Zj imax ^ dRi for some z = 1, . . . , N then VP(zj imax ) = 0, 
therefore 

(26) P(z hmax ) - P{x) > min | VP{y) | dist(z, -9P,) 

for all i = 1, . . . , N and a; G P«. Let Aj be a spherical cap of radius bdN^ 1 ^ 
contained in Ri. Since ||P|| < KdN~ l l d we obtain that 

(27) / dist(x, dRi)dfi d (x) > I dist(a;, dAi)d/i d (x) > Aj-^-, 
for some constant (3d- Thus, it follows from ( 1251) and ( 1261) that 



i N wnw 

£ P(w) > £ min | VP(y) | . 

j=l i=l 

Since by Theorem [2] we have ||"P|| < T]r d /(t + 1), using Lemma we arrive 
at 

1 * /• 

( 28 ) Tt J2 P (^,ma x ) > \\K\\P d (l-4dTl) / | VP(x)|^(l). 

^ 7=? ^ 
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Next we estimate the sum 

1 - 

-]T(PK max )-p(^)). 



N 

i=l 



Let 7 : [0, 1] — > Ri be a geodesic connecting Zi and -2j imax - We can write 

P{z i>max ) - P{ Zi ) = [ (VP( 7 (/i)),7'W) ^ > 0. 

By Lemma [3] (iii) the inequality ( , ~f'(h), VP(xj)) < holds true for all h 6 
(0, 1). Thus, we have 

/ 1 (VP( 7 (/ i )),7 / (^))^ < [ (VP(j(h)) - VP(x 4 ), -/(/»)) dfc 

< / I VP(7(/i)) - VP(xi)| | 7'(/i)| d/i < diamPi max | VP{x) - VP (a*) I- 
Using Lemma [2] we arrive at 

1 * f 

( 29 ) T7 5"l^(^max)-PN| < ||^|| 4^ / \VP(x)\d(l d (x). 

N l^t Js* 
Now we estimate the third sum in the left-hand side of (1241). Recall that 



Vi,e 



Hence, we obtain 



for \VP(Xi)\ > £ 

l[x uZi }{WP{xi)\/£) otherwise. 



1 N 1 

„ E( p (^) - p M) = n S ( p w - 

i=l i:\WP{xi)\<e 

Since Z{ and yi j£ are both in Ri, we can write an obvious estimate 
\P{yi,e) ~P{zi)\ < diam ifc max I VP(a;) I . 
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For each i = 1, . . . , N choose Wi G Ri such that | VP(iOj) | = max xgfi - | VP(x) | . 
Then 

1 £ |VPK)|< 

i:|VP(xi)|<e 

1 £ (|VPK)|-|VP(x 4 )|)+e+i £ (|VP(^)|-|VP(^)|). 

i:|VP«)|<£ i:|VP(xi)|>£ 

Thus, it follows from Lemma |2J that 



i:\SJ P{xi)\<e 

Hence, we arrive at 

N 



l - I VP K) \<s + 4:d V f | VP(x) I d Md (z) 



(30) 



^(Py-P(^)) < (e + 4d V J s jVP(x)\df, d (x)y 



It remains to estimate the sum 



1 N 

E(P(x,(P))-P( 2/i)£ )). 



The distance between Xj(P) and |/j i£ is less than 5 \\7l\\. Hence, 



1 N 

-^(P(x,(P))-P(^)) 

i=l 



5 \\1Z\\ v^->, i „ „/ x , 
<^^> max \VP(x) . 



Using again Lemma |2j we arrive at 



1 - /" 

( 31 ) ^E(P(x,)-P(^)) < 11^115(1 + 4^) / |VP(x)|d// d (x). 

Now, it follows from (124]), (1251). (125}. (I3"0l). and (|3"T]) that for Pe5fi 
1 - 

(32) ^^P(x,(P)) > \\n\\(l3 d (l-Adr ] )-Ad V ~(Adr ] + s)-5(l + Adr 1 )). 



i=l 



19 



Take 77 = p d /{24d),6 = f3 d /3 and e = d /12, where (3 d provided by (1271). 
Without loss of generality we may assume that pd < 1. Thus we get 

&(1 - Adr]) - Adr] - (4dr) + e) - 5(1 + Adrj) > 0, 

which together with ( 132]) imply f|22|) . 

It remains to show the separation property (123]) . Fix P G Q and i G 
l,...,iV. By the definition ( 12T1) . Xj(P) is in Thus to prove ( 1231) it is 
enough to show that dist(xj(P), dRi) > XdN^ 1 ^ for some constant A^. Recall 
that Xj(0) = X{, and Ri contains a spherical cap of radius r = b d N^ l / d with 
center at Xj. The main reason why Xi(P) is "far away" from the boundary 
dRi is because Xj(P) lies on the geodesic ■J[x i ,z l ], where G -R«, and 

dist(xi,Xj(P)) < (1 - 6)dist(xi,Zi). 

We will also use the fact that Ri is geodesically convex and contains a spher- 
ical cap of "big" radius with center at 2j. 

Denote by T x . the space of all vectors in M. d+1 orthogonal to Xi and let 
p : M. d+1 — y T x . be the orthogonal projection 

p(z) = Z - (Xi,z)Xi. 

As we have pointed out in Lemma [3] the image Si = p(Ri) is a convex subset 
in T Xi . Clearly, 

dist(x i (P),5P i ) > dist cuc (p(x i (P)),95 i ), 

where dist enc (z,dSi) stands for the Euclidean distance between point z and 
the set dSi in T x .. The point p(xj(P)) lies between the points p(xi) = and 
p(zi) on the line connecting them. Thus we have 

(33) \p{xi) -p{zi)\ = sm(dist(xi, Zi)), 
and 

(34) \p( Xi ) -p(xi(P))| < sm{{l-8)dist{x i ,z i )). 
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Moreover the fact that _R, contains the spherical cap A(xi,r) implies that 
(35) dist euc (p(x i ),9S'i) > sinr. 

Finally we note that the function dist euc (z,dSi) is concave on Si. Therefore 
by PD-flMD we get 

r 

&st(xi(P),dRi) > - sinr > \ d N~ l/d , 
which implies ( 1231) . Lemma H] is proved. 

□ 
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